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Introduction 


Volatile Removal Assembly (VRA) is a subsystem of the Closed Environment Life Support 
System (CELSS) installed in the International Space Station. It is used for removing 
contaminants (volatile organics) in the wastewater produced by the space station crews. The 
major contaminants are formic acid, ethanol, and propylene glycol. The VRA contains a slim 
packbed reactor (3.5 cm diameter and four 28 cm long tubes in series) to perform catalyst 
oxidation of wastewater at elevated pressure and temperature under microgravity conditions. In 
the reactor, the contaminants are burned with oxygen gas (CL) to form water and carbon dioxide 
(CCL) that dissolves in the water stream. Optimal design of the reactor requires a thorough 
understanding about how the reactor performs under microgravity conditions. The objective of 
this study was to develop a mathematical model to interpret experimental data obtained from 
nonnal and microgravity conditions, and to predict the performance of VRA reactor under 
microgravity conditions. 

Catalyst oxidation kinetics and the total oxygen-water contact area control the efficiency of 
catalyst oxidation for mass transfer, which depends on oxygen gas holdup and distribution in the 
reactor. The process involves bubbly flow in porous media with chemical reactions in 
microgravity environment. This presents a unique problem in fluid dynamics that has not been 
studied. Guo et al. (2004) developed a mathematical model that predicts oxygen holdup in the 
VRA reactor. No mathematical model has been found in the literature that can be used to predict 
the efficiency of catalyst oxidation under microgravity conditions. 

Mathematical Model 

A mathematical model has been developed in this study to simulate the oxidation efficiency of 
the VRA reactor. The longitudinal distribution of the oxygen mass flow rate in the reactor can be 
expressed as follows (derivation of the model is available from the author): 


M 02 — M 02-0 e 


r bP 0 2 v b 


( 1 ) 


where M Q2 is oxygen mass flow rate at the point x-distance from the reactor inlet, M Q2 0 is 
oxygen mass flow rate at the reactor inlet, r Q2 is rate of oxidation per unit contact area between 
oxygen and water, x V b is the average sphericity of oxygen bubbles (the ratio of bubble’s shortest 
diameter to longest diameter), r h is the average radius of bubbles, p m is the density of oxygen 
gas, and v h is the average velocity of bubbles. 


It is customary to use the term oxygen utilization ( U 02 ) in evaluation of VRA reactor 
performance. The expression for U 02 can be theoretically derived based on Eq. (1): 
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( 2 ) 


where L is the length of the reactor and R 02 -L is the remaining oxygen at the outlet of the reactor. 
The oxygen mass consumption potential in the reaction is then expressed as: 

02 = U 02 M O2 _ 0 t w (3) 
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where t w is the residence time of water in the reactor. The water residence time can be calculated 
based water injection rate, reactor size, packbed porosity, and oxygen holdup in the reactor. 
Based on the molar equivalency in the oxidation, the amount of organics i removed in the 
reaction can be expressed as: 


AM, 


A MpiWm, 

32 A, 


( 4 ) 


where W m is the molecular weight of organics i, and A, is molar number of oxygen utilization for 
each mole of organics i defined in the following oxidation reaction: 


tit [organics] + m ; 0 2 => p,- C0 2 + qi H 2 0, 


( 5 ) 


i.e., A = nij/rii. Typical values of A are shown in Table 1. 


The oxygen mass flow rate at the reactor inlet M O2 _ 0 is a design parameter. The density of 
oxygen gas ( p 01 ) can be calculated with gas law for ideal gas under given pressure and 
temperature. The sphericity of bubble (T^) may be affected by reactor geometry and gas holdup 
in the reactor. Due to lack of means of estimating its value, we used 0.5 throughout this study. 
The average flow velocity of oxygen bubbles ( v h ) can be estimated based on oxygen injection 
rate, reactor cross-sectional area, packbed porosity, and oxygen holdup. The oxygen holdup can 
be estimated with the mathematical model described by Guo et al. (2004). 


The rate of oxidation per unit water-oxygen contact area ( r Q2 ) depends on the total flux rate ( J ) 
at which volatile organics disperse (Taylor, 1953): 
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where D oi is the molecular diffusivity of organics i, F is called fonnation factor which is defined 

dC 

as tortuosity divided by porosity of packbed, — L is the concentration gradient of organics i in 

dy 

the direction perpendicular to the bubble surface, v is the average convection velocity, and r is 
the average pore radius. 

A mathematical model was developed in this study to estimate the size of gas bubbles in a 
flowing liquid (derivation of the model is available from the author): 
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iV WeC 


pL d n 




1.1548, 


g d m (.Pi - P g ) 


p L c d 


-v, 


( 7 ) 


where Nwec is the critical Weber number, p L and p are densities of liquid and gas, respectively, 
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v j is liquid velocity, and O is drag coefficient. If all other parameter values are known, the 
maximum bubble diameter d m can be solved numerically from Eq. (7). The critical Weber 

number was determined to be 979 based on matching the calculated bubble diameter to the data 
(air bubble size in water). By matching the calculated oxygen bubble diameter to the observed 
oxygen bubble length (4 to 6 inches) under 1-G conditions, we estimated that the drag coefficient 
in the VRA reactor to be about Cd = 3.9. 

Model Calibration 


The VRA reactor performance data obtained from ethanol oxidation an experimental test 
conducted in ground surface conditions were used to calibrate the model. The reactor influent 
contained 144 mg/L of ethanol. The reactor effluent contained 91 mg/L of acetic acid and no 
ethanol. Material balance calculations show that the oxygen utilization is 87%. This value was 
obtained from the model by assuming that the organics concentrations drop from the initial 
concentrations in the pore space to zero at the catalyst surface within a distance of a half pore 
radius, and adjusting convection velocity to 0.2957 times the longitudinal water velocity. 

Results 


Figure 1 presents calculated oxygen mass flow rate under 1-G and 0.0 1-G conditions. In this 
case, the oxygen utilizations under 1-G and 0.01-G conditions are 63.42% and 60.11%, 
respectively. The calculated oxygen utilization drops by 5.23% under micro gravity conditions. 

Conclusions 


A mathematical model has been developed in this study to estimate oxidation efficiency in the 

VRA reactor under microgravity conditions. The following conclusions are drawn. 

1. The remaining oxygen gas in the VRA reactor should decline exponentially with reactor 
length. The decline rate is directly proportional to rate of oxidation per unit area and bubble 
sphericity, and inversely proportional to bubble size, oxygen density and bubble velocity. 

2. Gravity affects oxygen utilization through changing several parameters including oxygen 
bubble size. More research efforts are required to further investigate the effects of 
microgravity on VRA reactor perfonnance. 

3. The currently VRA design is very conservative. The oxygen utilization is about 60%. 
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Table 1 - Baseline Concentration of Volatile Organics in the Waste Water 



Molecular 

Weight 

Liquid Molecular 
Diffusivity (cnr/s) 

Molar Number of 

Influent 

Organics Name 

Oxygen 

Utilization 

Concentration 

(mg/L) 

Acetaldehyde 

44.05 

0.00006953 

2.5 

2.74 

Acetone 

58.08 

0.00005924 

4 

4.41 

1 -butanol 

74.1 

0.00005924 

6 

1.5 

Ethanol 

46.07 

0.00006953 

3 

39.2 

Ethylene glycol 

62.07 

0.00005924 

2.5 

5.51 

Formaldehyde 

30.03 

0.00006953 

1 

5.89 

Formic Acid 

46.03 

0.00006953 

0.5 

49.9 

Methanol 

32.04 

0.00006953 

1.5 

7.94 

Propylene glycol 

76.1 

0.00005924 

4 

35.4 

1 -propanol 

60.1 

0.00005924 

4.5 

0.45 

2-propanol 

60.1 

0.00005924 

4.5 

1.6 

Urea 

60.06 

0.00005924 

1.5 

3.55 



Distance from the Entry Point (cm) 


Figure 1 - Model-Calculated oxygen mass flow rate under 1-G and 0.0 1-G conditions 
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